Kinetic growth of field-oriented chains in dipolar colloidal solutions 
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Experimental studies on the irreversible growth of field-induced chains of dipolar particles suggest 
an asymptotic power-law behavior of several relevant quantities. We introduce a Monte Carlo 
model of chain growth that explicitly incorporates the anisotropic diffusion characteristic of a rod- 
like object. Assuming a simple power-law form for the mean cluster size, S(t) ~ t z , the results of 
our model are in good agreement with the experimental measurements of the dynamic exponent z. 
Nevertheless, an alternative scenario, including logarithmic corrections to the standard power-law 
behavior, provides a better and more insightful interpretation of the anomalous dynamic exponent. 
In contrast to some experimental findings, we do not observe any dependence of the exponents on 
the volume fraction of particles <f>. Finite-size effects are also explored by simulating very long time 
evolutions or highly concentrated systems. Two different behaviors are found, namely, saturation 
and a crossover to a quasi one-dimensional regime. 

PACS number(s): 61.43.Hv, 82.70.Dd, 64.60. Cn 



I. INTRODUCTION 

The kinetic properties of the irreversible aggregation 
of particles has been the subject of great interest over 
the last decade In particular, several experimental 

and theoretical studies, as well as computer simulations 
have been devoted to understanding the behavior of col- 
loidal suspensions of dipolar particles, the so-called elec- 
trorheological and magnetorheological fluids [p|-|rT| . After 
applying an external electric or magnetic field, these sus- 
pensions experience a dramatic change in their rheologi- 
cal properties — i.e., a notable increase in the viscosity — 
which renders them particularly interesting systems from 
a technological point of view. 

The change in the rheological properties of dipolar sus- 
pensions upon the action of an external field, is mainly 
due to the aggregation of the colloidal particles, which 
form clusters of macroscopic size. These are usually 
linear chains — rods — oriented along the direction of the 
applied field, although for high enough concentrations 
of dipolar particles more complex structures may arise 
jl2],[l3|]. The overall spatial arrangement of the aggregates 
is very effective in hindering the fluid flow, conferring the 
suspension a solid-like texture. 

Multiple applications have been envisaged for these 
materials: lubricants, dampers, heat and light transmis- 
sion devices, etc. Nevertheless, their eventual manufac- 
ture still has to overcome different drawbacks, as for in- 
stance, the sedimentation of the clusters. At high fields, 
the aggregation is irreversible, the chains grow longer and 
never break as long as the field is present. Thus, the 
study of the kinetics of formation of the aggregates turns 
out to be an issue of practical importance. 

Due to its irreversible character, studies have most of- 
ten focused on the dynamic properties of the distribu- 



tion of clusters, following the approach developed to deal 
with irreversible cluster-cluster aggregation models (see 
Refs. and references therein). Within this frame- 

work, the quantitative features of the aggregation are 
commonly described in terms of the cluster-size distribu- 
tion as a function of time, n s (t), which is defined as the 
number of clusters of size s, per unit volume, present in 
the system at time t. Other relevant quantities are the 
density of clusters at time t, n(t) = X) s n s(i), anc ^ the 
mean cluster size S(t), defined by 



S(t) = 



(1.1) 



Experimental results and computer simulations of differ- 
ent cluster-cluster aggregation models |||| show that the 
asymptotic behavior at large times of the mean cluster 
size and the total number of clusters is a power law: 



S(t) ~ t z 



n(t) ~ t~ 



(1.2) 



The exponent z is the so-called dynamic exponent. In 
general, z depends on the dimension d of the space, as 
well as on the nature of the aggregation process. An- 
other relevant feature of the aggregation dynamics is 
that the characteristic quantities defined above are re- 
lated through the dynamic scaling hypothesis ]T^ |: 



n s (t) 



2 F[s/S(t)}, 



(1.3) 



where F(x) is a certain scaling function, independent of 
t and s. The validity of this relation has been systemat- 
ically verified in experiments and computer simulations. 

The dynamics of cluster-cluster aggregation can be 
theoretically described in terms of the Smoluchowski 
equation flloj ]. Smoluchowski proposed a kinetic equa- 
tion describing the temporal evolution of the cluster-size 
distribution n s (t): 



1 



dn s (t) 
dt 



= K(i,j)m(t)nj(t) - n s (t)J2K(s,i)m(t). 

i=l 

(1.4) 



i+j=s 



Here K(s, s') is the reaction kernel, giving the rate at 
which clusters of size s join clusters of size s' to form 
clusters of size s + s'. This equation is an example of 
a mean-field-type theory in which fluctuations are ne- 
glected. Moreover, its range of applicability is limited to 
low concentrations, where the assumption of exclusively 
binary collisions is valid jL5|. The results of this theory 
are expected to hold for dimensions higher than the up- 
per critical dimension (d c = 2 in this case fl6| ), above 
which fluctuations become irrelevant. Assuming that the 
reaction kernel is a homogeneous function of degree 7, 
that is, 



K{bs,bs') = b^K(s,s'), 



(1.5) 



it is possible to prove the relations Jl.2| ) and (1.3) ||, 
with a dynamic exponent given by 

(1.6) 



1 



' 7 



Most of the experimental and numerical studies of the 
anisotropic aggregation of rod-like clusters presented so 
far [[^Jq, |l0| , |l7| seem to indicate the same type of power- 
law for the asymptotic behavior of the mean cluster size 
S(t). For instance, Fraden et al. || reported a value 
of z = 0.60. Promislow et al. ||, on the other hand, 
measured a value of z ranging between 0.50 and 0.75 
for different values of the concentration and the dipo- 
lar interaction strength A (to be defined later on). The 
interpretation of these results relies on Smoluchowski's 
assumptions. 

In Smoluchowski's solution fll5|] , the reaction kernel is 
the product of an effective diameter — the collision cross 
section of two clusters — (R s + R s i), times an effective rel- 
ative diffusion coefficient (D(s) + D(s')), where R s and 
D(s) are the radius of influence and the diffusion coeffi- 
cient of a cluster of mass s, respectively [ fL8[ . By analogy 
with the Stokes-Einstcin relation for a single spherical 
particle of diameter a moving in a liquid, D cx a -1 fll9| , 
the effective diffusion coefficient is given by a power law 
of the cluster size, 



D(s) 



(1.7) 



Miyazima et al. |10[] proposed this type of power-law 
kernel for rod-like clusters of dipolar particles. One 
can estimate the effective cross section — the radius of 
influence — of a single magnetic particle, R\, by com- 
paring dipolar (Ud(r) oc m 2 /r 3 ) and thermal (fcsT) 
energies. If we introduce the dimensionless parameter 
A = m 2 /(a 3 fcsT), where m is the magnetic moment and 
a the diameter of a dipolar particle, the effective cross 
section can be expressed as R\/a ~ A 1 / 3 . The detailed 



form of the dipolar interaction is not expected to mod- 
ify the asymptotic behavior of the aggregation dynamics. 
Thus, roughly speaking, outside a spherical region of ra- 
dius R\ the relative motion is mainly diffusive, and only 
when one particle enters the sphere of influence of an- 
other they stick irreversibly. 

At low concentrations, it is reasonable to assume that 
rod-like clusters of dipolar particles essentially aggregate 
tip-to-tip. Miyazima et al. argued that, under such con- 
ditions, the cross section of a chain cannot depend on 
its total length ^TJ, and it is possible to approximate 
R s ~ Ri. The only source of dependence on the cluster 
size should thus come from the effective diffusion coeffi- 
cient, and the reaction kernel is approximately given by 
K(s,s') ~ s 7 + s' 7 . This form fulfills the homogeneity 
condition (h.5 [) an d yields therefore a dynamic exponent 
given by Eq. ( |1.6|) . 

Following this approach to the problem, the value of 7 
is the key parameter in order to interpret the experimen- 
tal data |^,^J|]. The dynamic exponent z is numerically 
computed from a log-log plot of the mean chain length 
S(t) as a function of time. The value obtained is then 



associated through Eq. (1.6) to a particular value of 7. 
In the absence of hydrodynamic interactions among the 
spheres in a chain, the effective mobility is inversely pro- 
portional to its length. Consequently, o ne w ould expect 7 
to be equal to —1 which, according to ( |1.6D , corresponds 
to z = 1/2. In general, this value of z deviates from the 
ones measured experimentally. A simple way of inter- 
preting this difference is to consider an ad hoc value of 
7 different from —1, which phenomenologically accounts 
for a more realistic mobility of a rod- like particle [pl| . 

In this paper we propose an alternative mechanism to 
explain those discrepancies in the measured value of z. 
The core of our proposal is a different view of the mobil- 
ity of a cluster of mass s. Fraden et al. || already pointed 
out that hydrodynamic interactions should increase the 
effective mobility of a cluster of mass s, making it larger 
than the mobility of a collection of s independent par- 
ticles. The latter approximation, stating D(s) ~ s , 
is strictly valid for clusters of hydrodynamically nonin- 
teracting spheres plfl . However, this is too strong an 
assumption in the case under consideration, where the 
magnetic spheres form relatively rigid chains. In this 
case, it is well-known that the hydrodynamic drag for 
the translational motion is anisotropic [ pl|j23| . The drag 
coefficients along the direction parallel to the axis of the 
chain and in the perpendicular directions, £|| and £j_ re- 
spectively, are approximately given by p3fl 



2irna 



ln(s) 



and 



£l = 2£y, 



(1.8) 



where n is the viscosity of the solvent and s is the num- 
ber of particles composing the cluster — its mass. The 
mobility of a cluster is defined as the inverse of the drag 
coefficient; hence, the constants Du and D±, characteriz- 
ing the diffusion parallel and perpendicular to the rods' 
axis, are given by |1| 



2 



Du 



ksT ln(s) 

h s 



and 



D 



D \\ 
~~2 



(1.9) 



That is to say, hydrodynamic interactions generate 
anisotropic diffusion coefficients, exhibiting logarithmic 
corrections to the previously considered power laws . In 
order to investigate the effect of the new diffusivities fll.Sj ) 
on the dynamics of the process, we have introduced them 
in a Monte-Carlo 24 computer model for the cluster- 
cluster aggregation of rod-like particles. 

By simply considering anisotropic diffusivities with 
logarithmic corrections, and assuming that S{t) ~ t z , 
we are able to recover a dynamic exponent z in good 
agreement with experimental results. However, a simple 
heuristic argument suggests a different functional form, 
namely, a power-law with logarithmic corrections. 

Let us assume that all clusters in the suspension have 
the same average length S. The average separation R 
between neighboring clusters in a suspension with ini- 
tial volume fraction <f> — the initial density of dipoles — can 
then be estimated to be 



1/d 



(1.10) 



In d = 1, only two neighboring clusters are able to ag- 
gregate. In general, this is the most likely event in low 
dimensions. In high dimensions, however, and due to the 
diffusive nature of their movement, any two clusters are 
equally likely to join, irrespective of their relative dis- 
tance, as they are quite invisible to each other. A mean- 
field type of behavior is thus expected, in the sense that 
only the density of clusters present n(t) is relevant, and 
not their spatial arrangement. 

As the movement of the clusters is essentially diffusive, 
the average displacement of a cluster of mean size S after 
the time interval t is given by 



(R 2 ) oc D{S) t, 



(1.11) 



where, as indicated in Eq. fll.9,) , D(S) ~ ]n(S)/S. Ag- 
gregation of two clusters will occur after a characteristic 
time T has been elapsed; a time interval long enough for 
the clusters to cover their relative separation. 

For a low-dimensional system, the characteristic time 
T is that required to cover the distance R separating near 
neighbors, i.e., 



Tex 



if 



1 g(2+d)/d 



D(S) <pi d ln(S) 



(1.12) 



On the other hand, in higher dimensions there is no 
such characteristic length scale and, consequently, there 
is no reason to expect that the previous expression holds. 
However, we can argue that a cluster browsing a volume 
of order unity will encounter n(t) clusters available to join 
to. The time needed for a cluster to cover such space 



will be thus proportional to nT, where, as above, T is 
the characteristic time of a single aggregation event. We 
have then 



Toe - 



1 1 



1 

oc — 



S 2 



n D(S) <f> ln(5) 



(1.13) 



where we have used that n ~ 4>/S. Note that both esti- 
mations yield the same result at d = 2, suggesting that 
the upper critical dimensio n of the pr oblem is d c — 2 [fl6| . 

The inspection of Eqs. (1.12) and ( 1.1 3| ) suggests that 
the functional dependence of the mean cluster size with 
time takes the form, for d < d r , 



S 



MS)] 



c 



W> 2/d ) c , 



where £ = d/(2 + d); and for d> d c , 
S 



MS)] 1 ' 3 



1/2 



(1.14) 



(1.15) 



Note that within this approach, we obtain logarithmic 
corrections to th e b ehavior reported in Ref. (To) for 
7 = — 1 (sec Eq. (|l.6|) ). These corrections to the usual 
asymptotic behavior in isotropic cluster-cluster aggre- 
gation could explain the anomalous dynamic exponent 
found for the aggregation of anisotropic rod-like clusters. 
The results of our numerical simulations exhibit a sur- 



prisingly good fit to the theoretical predictions (1.14) and 



(1.15), better than to the naive power law behavior. 

We have structured our paper as follows. In Section |H| 
we describe the technical details of our algorithm. Sec- 
tion [II deals with the properties of rod-like aggregation 



in two dimensions. In particular, we study the mean clus- 
ter size S(t) for a variety of values of the initial volume 
fraction <f>. For small to moderate values of <p and not very 
long execution times, we recover, assuming a pure power 
law form for S(t), values of z in agreement with some of 
the experimental findings. Importantly, in opposition to 
some claims in the literature ||, we do not observe any 
dependence of z on the volume fraction. On the other 
hand, we observe that the same sets of data can be fitted 
with higher accuracy to our predicted power law with 
logarithmic corrections. By allowing very high volume 
fractions or large execution times, we find that our model 
crosses over to two different behaviors, namely, satura- 
tion and a quasi one-dimensional regime. We believe that 
the discrepancies with the experimental results in Ref. Q 
might be due to a combination of this type of finite-size 
effects and logarithmic corrections. In Sec. IV we extend 
our model to three dimensions. We do not observe sig- 
nificant variations with respect to the two-dimensional 
aggregation. Finally, our conclusions are presented in 
Sec. M. 
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II. COMPUTER MODEL 

In our model we consider the irreversible aggregation 
of rigid rod-like clusters in H , with d = 2 and 3. Clus- 
ters are oriented along the Z axis. Simulations start at 
t = with a random distribution of Nq monomers (spher- 
ical particles of diameter a = 1) in a box of volume V 
with periodic boundary conditions. The initial volume 
fraction of monomers is defined by <j> = Nq/V. We have 
simulated a wide range of volume fractions <fi (0.001 to 
0.1), and some particular high values (up to 0.5). Clus- 
ters diffuse performing a free off-lattice random walk; 
that is, despite the dipolar interactions, we assume that 
the temperature in the system is high enough to provide 
a mainly diffusive character to the aggregation dynam- 
ics. Nevertheless, dipolar interactions are predominant 
at short distances. Therefore, when two clusters come 
close enough to each other, they stick irreversibly. 

In the Monte-Carlo algorithm, clusters are selected and 
moved a distance equal to one diameter a in a direction 
chosen randomly from a certain probability distribution. 
The movement is performed rigidly, preserving the ori- 
entation of the rod along the Z axis. When two clusters 
come within a distance "tip-to-tip" of one diameter (that 
is, when the distance between any two of their respective 
ends is less than or equal to a) they join, forming one sin- 
gle rod of mass — number of particles — equal to the sum 
of the masses of the colliding clusters. When two clusters 
approach "side-to-side" they repel each other. This is im- 
plemented by rejecting all possible movements leading to 
a "side-to-side" overlap of clusters. 

In fact, this procedure is equivalent to considering a 
radius of influence R s — a, and, according to the dis- 
cussion in the previous section, it mimics a dipolar in- 
teraction of strength A ~ 1. Different values of A could 
be in principle simulated by joining clusters irreversibly 
when they come within a distance ~ A 1 / 3 . Values of A 
between 1 and 30 (that is, radius of influence between 1 
and 3), within the range reported in experimental inves- 
tigations, provide analogous asymptotic results at large 
times, differing only in the transient regimes. 



p(s) ~ D c {s) 



The effect of the anisotropic diffusivities (L9) on the 
algorithm is taken into account in two steps: First, in 
the selection of the next cluster to be moved, and sec- 
ond, in the prescription for the direction of its tentative 
movement. In the first step, we choose a cluster among 
all present at a given time step, with a probability p(s ) 
proportional to their diffusivity. According to Eq. (1.9), 
we have 



P(s) 



ln(s) 



(2.1) 



However, given that the expression (l.S) is only suitable 
to describe the motion of a slender body (s 3> 1), the 
previous expression would be inappropriate for small s. 
Therefore, we actually select the cluster with a probabil- 
ity proportional to the corrected diffusivity 



ln(s )/s 
ln(s)/s 



for s < s 
for s > so 



(2.2) 



Here we are assuming that small clusters (s < sq) diffuse 
as if composed of hydrodynamically independent parti- 
cles. We have selected a cut-off mass s = 3; different 
values were also tested, yielding comparable results. 

We discuss now the method employed to sample the 
clusters. In cluster-cluster aggregation simulations, the 
sampling with a probability proportional to the mobility 
is usually performed according to the "rejection" algo- 
rithm ( |p|,p5| and Ref. [p6| , p. 151): A cluster is selected 
among all the present at a given time, with uniform prob- 
ability. Then, a random number 77, uniform in the inter- 
val [0, 1], is drawn. The cluster selected is accepted for 
movement if r\ < -D(s)/Z) max , where D(s) is the diffu- 
sion coefficient of the cluster considered, and D max is the 
maximum diffusivity of all the clusters present. Other- 
wise, the cluster is rejected, and both selection steps are 
repeated. In our simulations, however, we have chosen 
to implement a different sampling algorithm, the "alias" 
method for discrete distributions (see Ref. p3), p. 158, 
for a detailed description). Briefly, the algorithm works 
as follows: A cluster C, of size s, is selected uniformly 
among all clusters. Then, with a certain "aliasing" prob- 
ability p(s), the cluster is replaced by its alias C = F[C], 
of mass s' . The probability p(s) is chosen so that clusters 
with small Dc{s) (low probability of being selected) are 
frequently mapped to clusters with higher Dc(s'). The 
"alias" method is more costly in computer time; however, 
it provides a more accurate sampling of the probability 
density p{s). 

With the procedure described above, we implement the 
mass dependence of the mobility (1.9). However, we also 
have to take into account its anisotropy. Given Eq. ( [1.9] ), 
it is twice more likely for any rod to move along its axis 
that along any other perpendicular direction. We imple- 
ment this fact in a second step by selecting the direction 
of the trial movement from a probability distribution ful- 
filling this same anisotropy. Let us define the azimuthal 
angle 8 with respect to the axis of the rod. Then, in 
d = 2, the direction for the tentative movement of a 
chain of length s > sq is selected at random from the 
probability density 



P{8) 



V2 



1 



2tt 1 



(2.3) 



with < 8 < 2ir. In d — 3 we use instead the density 
1 \/2 1 



P(cp,8) = 



2tt ln(3 + 2V2) 1 + sin 



(2.4) 



with < <p < 2ir and < 8 < n. These distribu- 
tions ensure the necessary anisotropic condition P[8 = 
0, tt]/P[8 = 7r/2,37r/2] = 2, while being continuous in 
both angles, and easily simulated numerically. For clus- 
ters of mass s < sq we assume that the diffusion is 



4 



isotropic. Therefore, in this case we choose a random 
direction for the tentative movement from a uniform dis- 
tribution. 

The final implementation of our algorithm runs as fol- 
lows: (i) Each time step we select at random a clu ster of 
mass Si according to the prob ability density ( |2.2| ) and a 
direction Oi, according to ( |2 .3| ) [a pair (9i,ifi), according 
to (2.4), in three dimensions] for Sj > so; in case that 
Si < so, the direction is drawn from a uniform distribu- 
tion, (ii) The cluster is moved a distance a in the selected 
direction, and its position with respect to the neighbor- 
ing chains is analyzed, (iii) If it intersects "side-to-side" 
with another chain, the movement is rejected; otherwise, 
it is accepted. If the cluster intersects "tip-to-tip" with 
another chain, they join and form a single cluster, (iv) 
Finally, the time is incremented by an amount 



At 



N(t)D c ( Sl 



(2.5) 



where N(t) is the total number of clusters present at the 
corresponding time t, and Dc(si) is the effec tive diffu- 
sivity of the selected cluster, as given by Eq. (2.2). This 
choice of the time step effectively reproduces the real dy- 
namics of the system. 



III. AGGREGATION IN d — 2 

In our simulations in d = 2 we have mainly considered 
systems with volume ranging from 256x 128 to 1024x512, 
and initial number of particles between 524 and 67109. 
This corresponds to volume fractions ranging between 
4> ~ 0.001 and ~ 0.1. We have also considered the 
extreme case of high concentration (</> ~ 0.5) in order 
to investigate a possible crossover to a one-dimensional 
regime. In most cases, averages were performed over 100 
to 500 simulations. 
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FIG. 1. Aggregation with cluster diffusion that is isotropic 
and proportional to s _1 . Initial concentration (f> = 0.019. a) 
Mean cluster size as a function of time, b) Finite-size scaling 
of the cluster size distribution as function of time, for fixed 
lengths s = 10, 20, 30, 40, 50, 60, 70, and 80. 



for a value of z = 0.50, which is again in agreement with 



(1.6) 



A. Isotropic diffusion 



B. Anisotropic diffusion: dilute regime 



Firstly, in order to check our algorithm, we have simu- 
lated the aggregation according to the prescription given 
in Ref. pp| , i.e., the diffusivity of a rod is isotropic and 
inversely proportional to its mass, D(s) ~ s _1 . 

In Fig. 0(a), we represent the mean cluster size S(t) 
obtained from simulations of a system of size 1024 x 512, 
with an initial number of particles iVo = 10000. These 
values correspond to an initial volume fraction <f> = 0.019. 
We observed a clear power law regime covering more than 
three decades in time. The least-squares fitting of the 
curve yields an exponent z = 0.50. As expected, this 
result matches the dynamic exponent predicted by fll.6| ) 
for 7 = — 1. In Fig. jj(b) we have tested the finite-size 
scaling relationship ( |F3| ) by plotting s 2 n s (t) as a func- 
tion of the rescaled time t/s 1 ^ 2 for different values of the 
cluster size s. The best collapse of the plots is obtained 



We now discuss the results obtained by implementing 
into our algorith m th e anis otro pic diffusion prescribed in 
d = 2 by Eqs. and (|]|). First of all, motivated 

by the results in Ref. || , we have duplicated the param- 
eters in that experiment, i.e., we have chosen a system 
of volume 1024 x 512 and an initial number of parti- 
cles Nq = 5000. These conditions yield a volume frac- 
tion <f> ~ 0.01, as in the experiments. Figure ||(a) shows 
the mean cluster size as a function of time, in a double- 
logarithmic plot. The behavior of S(t) at late times can 
be fitted to a power law extending close to three orders of 
magnitude. A least-squares fitting provides an exponent 
2 = 0.61, in complete agreement with the experimental 
findings in ||. Similarly, in Fig. ||(b) we represent the 
finite-size scaling hypothesis ( |1.3[ ). The good collapse 
shown in this plot corroborates the value of the dynamic 
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FIG. 2. Aggregation with cluster diffusion that is 
anisotropic and proportional to ln(s)/s, as defined in 
Eqs. (2.2) and (2.3). Initial concentration (f> = 0.009, as in 
Ref. pi. (a) Mean cluster size as a function of time, (b) Fi- 
nite-size scaling of the cluster size distribution as function of 
time, for fixed lengths s = 10, 20, 30, 40, 50, 60, 70, and 80. 



exponent z — 0.61. 

The evidence exposed in the previous figures seems to 
lend support to the hypothesis of a simple power law be- 
havior. However, the situation is not completely clear, 
as we show in F igure ||. Fig. ||(a) tests the theoretical 
prediction (1.14). Plotting S(t) as a function of tln5 
in logarithmic scale we obtain a slope £ = 0.51, very 
close to the expected value 1/2 in d = 2. Moreover, 
we remark that the goodness-of-fit in this case, as mea- 
sured by the Pearson's r coefficient p7[ , is higher than 
the one obtained from the linear regression in Fig. ||(a). 
Indeed, least-squares fittings covering the last three or- 
ders of magnitude in the abscissae of Figs. ||(a) and ||(a), 
yield the values r = 0.99960 and r = 0.999948, respec- 
tively. In Fig. ||(b) we check the finite size scaling hy- 
pothesis by plotting s 2 n s (t) as a function of the rescaled 
time tlns/s 1 ^. The best collapse again results from a 
value of C = 0.51. 

In the simulations presented above we have been able 
to match the observed value of z by selecting the partic- 
ular value of <fr reported in the experiments. It has been 
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FIG. 3. Same data as in the previous figure, (a) Mean 
cluster size as a function of t In S(t). (b) Finite-size scaling of 
the cluster size distribution. 



argued, however, that the value of z may depend on the 
volume fraction tf> M . In order to verify the accuracy of 
this statement, we have performed simulations for several 
values of tj>, ranging from very dilute ((f) ~ 0.001) to mod- 
erate concentrations (4> ~ 0.1), in a system of fixed size 
1024 x 512. Figure ^(a) plots the mean cluster size S(t) 
as a function of t for the different volume fractions con- 
sidered. The slope of the different graphs depicted in this 
figure does not seem to depend on </>. This fact becomes 
even more clear in Fig. |](b), where we have plotted S(t) 
as a function of the rescaled time <pt. We observe that 
in this case all the plots collapse onto a universal func- 
tion, independent of the volume fraction. The collapse 
of the different graphs is also shown in Fig. |5|, now as a 
f unct ion of the rescaled quantity (f>t\nS. Note that Eq. 
( 1.14 ) leads directly to the required scaling factor of <f> in 
d = 2. As in Figure |^, the collapse is statistically better 
in terms of this new rescaling than for the single power 
law. 

The collapsed plots allow us to select a common scal- 
ing region for all of them, from which we can extract 
slopes for the different values of cj> that are directly com- 
parable. An average slope can thus be defined. In the 
case of Fig. ||(b) (single power law interpretation), indi- 



6 



2.5 
2.0 
1.5 

-t-o 

CO 

o 1.0 

bC 
O 

~ 0.5 



0.0 



a) 


4> = 0.064 


1 1 

////// ~ 






<y \ : 

= 0.001 







0.0 



2.0 



4.0 



6.0 



logio t 





2.5 r 




2.0 - 




1.5 - 




CO 




o 


1.0 - 


bJO 








0.5 - 




0.0 - 




-1.0 1.0 

logio(^) 

FIG. 4. Mean cluster size for anisotropic aggregation, for 
several values of the volume fraction. From top to bottom, 
4> = 0.064,0.032,0.016,0.008,0.004,0.002, and 0.001. (a) As 
a function of time, (b) As a function of the rescaled time <j>t. 
The average slope obtained from the common scaling region 
is 2 = 0.61. 



vidual slopes range between 0.59 and 0.63. From them, 
we obtain an average exponent z = 0.61, in accordance 
with our previous estimate and the experimental results 
in Ref. Q. In addition, from Fig. || (power law with log- 
arithmic corrections), we obtain slopes between 0.49 and 
0.53, yielding an average value of ( = 0.51. 



C. Anisotropic diffusion: saturated regime 

In view of the results presented in the last section, 
we naturally conclude that the slope does not depend 
on the volume fraction and that the disagreement with 
prior experimental findings Q may be a consequence of 
both logarithmic corrections and finite-size effects. In ex- 
periments or computer simulations, finite-size effects can 
cause the system to crossover, at large times, to a satu- 
rated regime in which any predicted scaling behavior is 
lost. We have investigated this issue by simulating sys- 
tems which were allowed to evolve for very long times — 
up to ten times longer than in our previous simulations. 
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FIG. 5. Same data as in the previous figure, now as a func- 
tion of the rescaled quantity <j>tlnS. The average slope ob- 
tained from the common scaling region is £ = 0.51. 




2.0 4.0 

log 10 [tln£(i)] 

FIG. 6. Mean cluster size as a function of tin S(t), for 
anisotropic aggregation in the saturated regime. The plots 
correspond to a fixed volume fraction <j> — 0.019 and differ- 
ent system sizes. The onset of the plateau and its height are 
functions of the system volume. 

Figure || depicts the mean cluster size computed for dif- 
ferent system sizes, while keeping the same fixed initial 
volume fraction, = 0.019. In Fig. ^, we plot S(t) for 
a fixed value of the system size 512 x 256, and different 
values of cf>. The chief feature of these plots is the onset 
of a plateau whose location and height appears to be a 
function of the system volume V. This flatter region is in- 
dicative of a considerable slowing down in the dynamics. 
Finite size effects, unavoidable at relatively large times, 
can corrupt the interpretation of any expected scaling 
behavior. 



D. Anisotropic diffusion: quasi one-dimensional 
regime 

We have established the existence of deviations from 
the dilute asymptotic behavior due to saturation when 
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FIG. 7. Mean cluster size as a function of the rescaled 
quantity <jrt]nS(t), for anisotropic aggregation in the 
saturated regime. The plots correspond to a fixed 
system size 512 x 256 and different volume fractions 
4> = 0.0048, 0.0095, 0.0191, and 0.038. The collapse of the 
curves shows that both the onset and the height of the satu- 
rated plateau depend only on the system size. 



the system evolves up to very long times. In addition, 
in a highly concentrated system, one can also observe 
a crossover to a quasi one- dimensional regime, as first 
remarked by Miyazima et al. |10| . At the late stages 
of evolution in a very concentrated system, the largest 
chains attain a length which is of the order of the system 
size. Their effect is therefore to divide the plane into 
one-dimensional strips. The aggregation of the remaining 
smaller clusters is thus mainly restricted to occur within 
these strips [gj, becoming effectively a one-dimensional 
process. 

In a truly on e-dim ensional model of aggregation, ac- 
cording to Eq. (1.14), we expect the mean cluster size to 
satisfy S(t) ~ [<p' z t lnS(t)} 1 / 3 . For the sake of complete- 
ness, we have adapted our model to simulate aggregation 
within a line. The only relevant change in the algorithm 
concerns the way in which the direction of the tentative 
movements is selected — either to the right or to the left, 
with equal probability. The results of the simulations are 
plotted in Fig. |[ This figure corroborates both the cor- 
rect scaling of time with 2 (in contrast to (f> for d > 2), 
and an exponent £ = 0.33. 

We next present results from simulations in d = 2 of a 
system exhibiting quasi one-dimensional behavior. The 
data in Figure ^ correspond to a volume fraction <j) = 0.48 
in a system of volume 1024 x 512 pq| . After a transient 
regime, we observe that the slope stabilizes in a value 
£ ~ 0.32, clearly smaller than the one found for less con- 
centrated systems, and in excellent agreement with the 
exponent ( — 1/3 expected in d — 1. We note that both 
the slope and the crossover time do not depend on the 
volume of the system for a fixed value of 0. The results 
obtained from simulations in a volume 512 x 256 (not 
shown) confirm this statement. 

We have further investigated the difference between the 
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FIG. 8. Mean cluster size as a function of the rescaled 
quantity 4> 2 tlnS(t), for anisotropic aggregation in d = 1. The 
plots correspond to a fixed system size 10000 and different 
volume fractions (f> = 0.02, 0.03, and 0.04. The average slope 
obtained from the common scaling region is £ = 0.33. 
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FIG. 9. Mean cluster size as a function of tin S(t), for 
anisotropic aggregation in d — 2, in the quasi one-dimensional 
regime. The dashed lines have slopes ( — 0.32, at larger times, 
and 0.50 at intermediate times, corresponding to the asymp- 
totic dilute regime. 



saturated and the quasi one-dimensional regimes. An 
effective way of doing so is to consider the movement 
of the individual clusters. As discussed above, in the 
quasi one-dimensional regime, average-sized clusters are 
confined to move within the narrow strips delimited by 
the largest chains. Therefore, it is more likely for a cluster 
to move along the direction parallel to the if-axis, than 
in any other direction perpendicular to it. This effect 
can be quantitatively assessed by estimating the relative 
frequency with which clusters move in a given direction 9. 
To this end, we define the jump orientation density F(6) 
as follows: At the last stages of the simulations (usually 
the last 25% of the time alloted for the run) we keep 
track of the direction 9 of all the accepted movements 
performed by the clusters. F(9)d6 is then defined as 
the probability that any of those movements is directed 
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FIG. 10. Jump orientation density F(9) for two systems of 
size 512 x 256, at two different volume fractions. The dilute 
system (<j> = 0.0048) matches the a priori angle distribution, 
Eq. ( pTs) ) (full line). The concentrated system {<j> = 0.38), in 
the quasi one-dimensional regime, shows two sharp peaks at 
6 = and 6 = tt. 



along a direction included in the interval [9,9 + d0\. At 
late times in the evolution of the system, most of the 
clusters have a length larger that the cut-off sq. If there 
is no bias in the direction taken by the accepted steps, 
we expect that the fu ncti on F(9) will match the a priori 
angular distribution (2.3). In the quasi one-dimensional 
regime, on the other hand, given that the clusters move 
with higher probability along the Z-axis, we expect F(9) 
to show anomalous maxima at 9 — and 9 = n. 

In Figure [l(] we plot F(9) as computed from two sys- 
tems of size 512 x 256, with very different values of <f>. 
The dilute system (<f> = 0.0048) is well inside the sat- 
uration regimen. As we can check from this plot, its 
angular distribution of jumps matches quite closely the 
a priory distribution of directions, plotted as a reference 
in full line. In the case of the highly concentrated sys- 
tem ((j) = 0.38), we can observe the two peaks in F(0), 
characteristic of the quasi one-dimensional regime. The 
secondary peaks at 9 = ir/2 and 9 — 3n/2 are a spurious 
result of the presence, even at the large time considered, 
of some clusters of size s < sq with an isotropic diffusiv- 
ity If they are removed from the statistics, the secondary 
peaks disappear. 



IV. AGGREGATION IN d = 3 

Our two-dimensional model can be easily extended 
to d = 3. In this case, the random directions o f th e 
jumps are sampled according to the distribution (2.4). 
Simulations have been performed for systems of size 
512 x 64 x 64. Results correspond to averages over 100 
realizations. 

In Fig. [ll] we plot our results for the mean cluster 
size S(t), for systems with an initial volume fraction 
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FIG. 11. Mean cluster size as a function of the rescaled 
quantity (f>tlnS(t), for anisotropic aggregation in d = 3. The 
volume fractions plotted are (f> = 0.0048, 0.024 and 0.048. 



(f = 0.0048, 0.024, 0.048. The average slope in the scaling 
region yields an exponent equal to the one obtained for 
two-dimensional aggregation, namely £ = 0.51. This fact 
is, however, not surprising: Since we are above the crit- 
ical dimension of the problem, d c = 2, we expect to find 
exponents independent of the dimensionality. As shown 
in the figure, the collapse of the mean-cluster size for dif- 
ferent volume fractions as a function of the rescaled time 
(f>t\nS(t), also holds in d = 3. 



V. CONCLUSIONS 

In this paper we have investigated the dynamics of 
the irreversible aggregation of linear rigid chains oriented 
along a preferred direction. A possible physical realiza- 
tion of such system would be the aggregation of dipo- 
lar particles in presence of a strong external field, ori- 
ented along the Z axis. We have proposed a Monte-Carlo 
model, whose key ingredient is a corrected anisotropic 
diffusivity of the clusters, as expressed in Eq. (1.8). This 
diffusivity exhibits a characteristic logarithmic correc- 
tion due to hydrodynamic interactions. A simple heuris- 
tic argument suggests that these logarithmic corrections 
should also emerge at a macroscopic scale, reflected in 
the asymptotic behavior of various quantities such as the 
mean cluster size. 

We have performed extensive simulations in d = 2, in 
a variety of system sizes and initial volume fractions. As- 
suming a simple power law behavior for the mean cluster 
size, S(t) ~ t z , our results yield a dynamic exponent in 
the scaling region z ~ 0.61. This value is in excellent 
agreement with previous experimental works || . For di- 
lute systems, the value of the dynamic exponent seems 
to be independent of the initial volume fraction. 

Alternatively, a better fit of our data is obtained in 
terms of the functional relationship S (t) / [\n(S (t))]^ ~ 
t'' , which explicitly incorporates logarithmic correc- 
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tions. Within this approach, the anomalous value of 
the dynamic exponent characterizing the aggregation of 
anisotropic rod-like clusters is, in fact, the consequence 
of the logarithmic corrections superposed to the usual 
power law behavior of isotropic cluster-cluster aggrega- 
tion. 

Long-time evolutions for any value of the volume frac- 
tion <f> lead to a saturated regime, characterized by a 
drastic slowing-down of the dynamics. The mean clus- 
ter size shows a flat region, whose height and onset are 
exclusively functions of the system volume. 

In highly concentrated systems, we observe a crossover 
to a quasi one-dimensional regime. Characteristic of this 
transition is a sharp change in the slope of the mean 
cluster size, yielding a value ( = 0.32. To support these 
findings, we have carried out simulations in d = 1 , where 
the mean cluster size exhibits a similar behavior with ( — 
1/3. In the quasi one-dimensional regime, clusters have 
a higher tendency to move along the direction parallel to 
their axis, than in any direction perpendicular to it. This 
point is made clear by analyzing the jump orientation 
density F{9). 

The extension of our model to d — 3 shows no sig- 
nificant variations with respect to the case d = 2. This 
result is to be expected, given that the critical dimension 
of the problem appears to be d c = 2. 

To conclude, we would like to remark the potential ap- 
plications of our model, especially the possibility of char- 
acterizing the aggregation of dipolar colloidal suspensions 
in different geometries of technological importance. For 
example, in a porous medium, due to the effective re- 
duced dimensionality, we would expect a slowing down of 
the aggregation rate, resemblant to the aforementioned 
quasi one-dimensional regime. 
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